Quantum simulation of many-body spin interactions with ultracold polar molecules 
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We present an architecture for the quantum simulation of many-body spin interactions based 
on ultracold polar molecules trapped in optical lattices. Our approach employs digital quantum 
simulation, i.e., the dynamics of the simulated system is reproduced by the quantum simulator 
in a stroboscopic pattern, and allows to simulate both coherent and dissipative dynamics. We 
discuss the realization of Kitaev's toric code Hamiltonian, a paradigmatic model involving four- 
body interactions, and we analyze the requirements for an experimental implementation. 
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I. INTRODUCTION 

The realization of quantum simulators - devices which 
can replicate the dynamics of other quantum systems 
[l|, Q - is currently one of the most exciting topics in the 
field of ultracold quantum gases [3]. This is particularly 
relevant in areas where classical simulation methods have 
proven to be inadequate due to the exponential growing 
Hilbert space dimension, such as in frustrated quantum 
magnets, where the system sizes that can be studied us- 
ing exact diagonalization methods are typically limited 
to less than 50 spins [4]. 

Ultracold polar molecules are particularly promising 
candidates for the quantum simulation of spin models due 
to their long coherence times and strong electric dipole 
interactions at distances compatible with optical address- 
ing. In the past, several theoretical proposals have been 
made for the realization of quantum magnetism in these 
systems (EI-IToj. However, these proposals are challeng- 
ing to extend to higher order many-body interactions as 
these terms arise within a perturbation series and thus 
become exponentially weaker the more particles partici- 
pate in the interaction [El, Q. On the other hand, such 
spin models with many-body interactions have recently 
received great attention in the context of Kitaev's toric 
code Hamiltonian which has interesting topological 
properties, and for the generation of cluster states, which 
are relevant for measurement-based quantum computing 

In this article, we describe an architecture for an effi- 
cient non-perturbative simulation of many-body spin in- 
teractions. Our approach is related to a recent quantum 
simulation proposal based on strongly interacting Ryd- 
berg atoms [13, 14]; however, due to the inherent limita- 
tions to Rydberg excitations due to their fast radiative 
decay, an implementation with ultracold polar molecules 
within their electronic ground state manifold might of- 
fer some beneficial aspects. Our approach to quantum 
simulation will be a digital one, i.e., the dynamics of the 
simulated system arises stroboscopically at discrete time 



intervals In particular, the effective dynamics will be 
created by applying sequences of microwave pulses cou- 
pling rotational excitations of the polar molecules, realiz- 
ing sequences of single-body and two-body quantum logic 
gates. We first describe the setup composed of ultracold 
polar molecules trapped in optical lattices, followed by a 
discussion on the implementation of a single many-body 
interaction. Besides coherent time evolution, the quan- 
tum simulator is also capable to incorporate dissipative 
dynamics, including efficient preparation of the ground 
state. We will then generalize these concepts to the full 
lattice system and discuss the sources and consequences 
of residual imperfections. Finally, we describe the exper- 
imental requirements to implement our quantum simula- 
tion architecture. 



II. SETUP OF THE SYSTEM 

We consider a two-dimensional system of ultracold po- 
lar molecules in their rovibrational ground state [HI, Il6j 
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FIG. 1: Proposed experimental setup. Polar molecules are ar- 
ranged in the intensity minima of a two-dimensional optical 
lattice with one molecule per lattice site. Molecules on neigh- 
boring lattice sites take different roles as control or ensemble 
molecules. 
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and loaded into an optical lattice [l^, see Fig. [H We fo- 
cus on the case where the lattice potential is deep enough 
to suppress any tunnehng between different sites. Fur- 
thermore, we assume the molecules to be initially pre- 
pared in a well-defined hyperfine state [18]. Here, we are 
also interested in a setup where a single control molecule 
is surrounded by four ensemble molecules^ in a way that 
the former can be manipulated independently without 
affecting the other ones. Such a controllability is most 
readily achieved if the molecules can be addressed indi- 
vidually using optical fields [l9|. Note that the distinc- 
tion between control and ensemble molecules is a purely 
logical one, i.e., both types can be realized using a single 
molecular species. 

The relevant level structure of the molecules is shown 
in Fig. [21 where we focus on the lowest four rotationally 
excited states. We consider the driving of rotational tran- 
sitions by microwave fields. By combining linearly and 
circularly polarized microwaves together with the electric 
quadrupole interaction coupling the nuclear spin to the 
rotation (S) [20], it is possible to selectively drive transi- 
tions between hyperfine states [l8|. 

Additionally, we are interested in the situation where 
a strong linearly polarized microwave field is reso- 
nant with the transition between the states | J = 2) and 
I J = 3). The quantization axis is defined by the po- 
larization, so the dynamics will be constrained to states 
with vanishing projection of the angular momentum on 
the quantization axis, i.e., Jz = 0. In the dressed frame 
of this driving, the state |-) = (| J = 2) - | J = 2))/v^ 
will pick up an effective permanent dipole moment of 
de = 3(i/\/70, where d is the bare electric dipole mo- 
ment of the molecule [lo[. Additionally, we consider a 
weak microwave driving Q of the two-photon transition 
between the state |0) and the dressed state | — ). 

Finally, the position of the |— ) manifold can be shifted 
relative to the | J = 0) manifold using optical fields result- 
ing in a differential ac Stark shift [21| , which can shift the 
molecules out of resonance of the microwave field Q. By 
confining the optical potentials to individual lattice sites, 
single molecule addressing on a lengthscale of 500 nm can 
be realized [l9| . Further improvements might be achieved 
using sub- wavelength addressing techniques using strong 
field gradients j22| or electromagnetically induced trans- 
parency (23| . 

Essentially, the setup presented here allows to selec- 
tively create rotational excitations, depending of the po- 
sition of the molecule and its hyperfine state. Conse- 
quently, the dynamics is effectively constrained to three 
states: two nuclear states corresponding to the | J = 0) 
manifold (named |0) and |1) hereafter), and one state in 
the 1+) manifold (|2)). The benefit of this setup based 
on microwave dressing compared to ones involving static 
electric fields lies in the strong suppression of dipolar flip- 
flop terms coupling different J states, thus eliminating a 
potential error source for the quantum simulator. 
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FIG. 2: Relevant part of the internal level structure. Depend- 
ing on the polarization of microwave fields (7r,cr^), hyperfine- 
preserving or hyperfine changing transitions can be driven. 
A strong microwave field couples the J = 2 manifold to 
the J = 3 manifold, resulting in dressed states that can be 
accessed by an additional two-photon microwave driving Q. 
For the purpose of the quantum simulator, only the states |0), 
and 1 2) are important. 



III. SIMULATION OF THE TORIC CODE 

The specific model we want to outline a quantum sim- 
ulator for is Kitaev's toric code [11]. It serves as paradig- 
matic model of a large class of so-called stabilizer Hamil- 
tonians [i^, [lH, whose ground states can be found by 
local minimization of the energy. The toric code Hamil- 
tonian is given by 



(1) 



with the "plaquette" 



operators Ap 

(i) 



n 



and the 



"vertex" operators By = Yliev containing Pauli oper- 
ators representing four-body spin interactions, see Fig.O 
On a torus, the ground state of the toric code is four- 
fold degenerate and corresponds to the +1 eigenvalues 
of all mutually commuting operators Ap and By and ex- 
hibits topological order. Plaquettes or vertices with a — 1 
eigenvalue are quasiparticles called magnetic charges or 
electric charges^ respectively, and always occur in pairs 
with a string operator connecting them Due to the 
non-commutativity of ax and a^, moving a quasiparticle 
of one type around a second one of the other type will re- 
sult in a minus sign applied to the quantum state; hence, 
the quasiparticles of the toric code are abelian anyons. 

These exotic properties have led to many proposed re- 
alizations within ultracold quantum gases or condensed 
matter systems [El, HH, [3, l26l-[29l ^ bu t despite experimen- 
tal success in small systems [30|, l31| , a large-scale many- 
body implementation revealing the topological proper- 
ties is still lacking [32|. In the following, we will discuss 
the implementation of the toric code within our quantum 
simulator architecture. As all terms of the Hamiltonian 
are mutually commuting, the dynamics of the toric code 
can be constructed by sequential implementation of the 
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stronger than the dipole-dipole interaction, 



Vdd 



J_4 

47r£o CL^ ' 



(3) 



where a is the separation between the control and the 
ensemble molecule. After the pulse, the molecules ex- 
perience the dipole-dipole interaction for the time t^^ = 
Trh/Vdd^ following a second microwave pulse transfering 
the population back to the rovibrational ground state. 
Then, up to local rotations which can be canceled by 
appropriate additional driving fields, the system has re- 
alized a conditional phase gate of the form 



(4) 



FIG. 3: Two-dimensional lattice setup for the toric code in- 
volving control (red) and ensemble molecules (blue) . The en- 
semble molecules taking part in plaquette operators Ap and 
vertex operators By are colored accordingly. 



where c and e refer to the control and ensemble molecule, 
respectively. 



2. Digital simulation procedure 



plaquette terms and vertex terms, i.e., 

U = ex.p{-iHt) = Y[^^p{^EoApt/h)Y[^^p{iEoByt/h). 

p V 

(2) 

The form of the plaquette interactions Ap and the ver- 
tex interactions By is identical up to a global rotation 
interchanging ax and a^; hence, we focus first on the 
implementation of a single By term, with all remaining 
terms of the Hamiltonian to be realized in an analogous 
way. 



A. Single plaquette interaction 

As a crucial step, we now make use of the control 
molecules to mediate the vertex interaction. During this 
process, we require the implementation of single qubit 
rotation between the logical states |0) and Such 
single qubit gates between nuclear spin states have al- 
ready been experimentally realized based on the electric 
quadrupole moment of the nucleus [l8|. Furthermore, 
our architecture will involve a two-qubit controlled phase 
gate between two molecules, which we will present in the 
following. 



1. Controlled phase gate 

The two logical states |0) and |1) do not exhibit any 
dipole-dipole interactions, hence, it is necessary to first 
transfer the population from the |1) state to the |2) state 
by a two-photon microwave pulse. Doing this simulta- 
neously in the control molecule and a single ensemble 
molecule requires the microwave driving Q to be much 



We will now turn to the quantum simulation of the 
vertex term By in the toric code Hamiltonian ([1]). Dur- 
ing a single timestep r of our digital quantum simula- 
tor, a single vertex undergoes the time evolution U = 
ex.-p{iEoByr). Such a dynamics is realized by a three- 
fold sequence, provided the control molecule is initially 
in |0)^'^'^ [m, (i) A sequence of quantum gates maps 
the eigenvalue ±1 of the operator By acting on the en- 
semble molecules is mapped onto the states |0) and |1) 
of the control molecule, (ii) A single qubit rotation of 
the form Uz{(t>) = exp(i(/)<ji^^) is applied to the control 
molecule, (iii) The mapping of step (i) is undone by the 
inverse gate sequence. The total gate sequence for the 
simulation of the time evolution reads 

U = exp{i(t)By) (5) 

e e 

where U^^^^ — exp(z<j^7r/4) is a 7r/2 rotation of the 
control molecule. This sequence involves a series of 
four controlled-phase gates Ucp acting on the control 
molecule and each ensemble molecule sequentially, which 
effectively represents a multi-qubit gate, where the con- 
trol molecule conditionally manipulates the ensemble 
molecules, sometimes denoted as a Controlled-Phase^ 
quantum gate. This sequential operation offers maxi- 
mum speed of the gate sequence as all operations can be 
carried out resonantly. On the other hand, conceptually 
slightly simpler but slower gate sequences based on direct 
multi-qubit gates [33| also exist [13, 14]. 

The effective energy scale Eq of the simulated Hamil- 
tonian depends on the phase (f) written onto the con- 
trol molecule during each timestep r, i.e., Eq = h(j)/r. 
While for the toric code involving only mutually com- 
muting operators, cj) can be arbitrarily large, models with 
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non-commuting degrees of freedom mandate the intro- 
duction of a Trotter expansion of the form exj[){—iHT) = 
Yl-ex.p{—iHir) + O(r^), requiring <C 1 to reproduce 
the desired dynamics with high accuracy. 

3. Dissipative state preparation 

As aheady mentioned, the ground state of the toric 
code Hamiltonian ([1]) can be found by locahy minimiz- 
ing the energy of the stabihzer operators Ap and By. 
Therefore, it is possible to engineer the dynamics such 
that energy is constantly removed from the system and 
the ground state arises as the final state of the dynam- 
ics As a crucial element, the dissipation 
of energy requires the presence of an incoherent process 
in the dynamics; here, we are interested in the situation 
where the control molecule can be incoherently pumped 
from the state |1) to the state |0). However, such a dis- 
sipative step is not as straightforward to realize as with 
atoms since the radiative decay of molecules to many vi- 
brational states makes optical pumping very challenging, 
although there has been considerable progress in laser 
coohng of molecules (sH, • Alternative routes to the 
desired dissipation process can be realized b y c oupling 
molecules to microwave stripline resonators [38|, [39| or 
by hyperfine-preserving STIRAP processes to dissociate 
molecules back into atoms, followed by optical pumping 
of the atoms and reversing the STIRAP process. 

For the implementation of dissipative state prepara- 
tion, we employ the gate sequence 

U = UcNOT,^Ul%UcpUl%, (6) 

where UcNOT,i is a controlled-NOT gate with the control 
molecule as a control qubit and the ith ensemble molecule 
acting as a target qubit. This two-body gate can be easily 
constructed from the controlled phase gate Ucp using 
additional single qubit rotations. The action of the gate 
sequence U can be understood by looking at the states 
|±,A) having the property (By) = ±1 (here, A encodes 
the information about the state that is independent of 
By). The chosen gate sequence will yield 

C/|0)J±,A), = -L(|o)J+i, A), + |1)J+1, A)J. (7) 

Performing the dissipative step |0)^ results in the 

state 1+1, A) for the ensemble molecules, independent of 
A. In general, the dynamics created by U and subsequent 
dissipation can be described in terms of a discrete time 
quantum master equation in Lindblad form, 

p(t + r)=cp(t)ct-i{ctc,p(t)}, (8) 

with the quantum jump operator c = iax\l — By)/2. 
This jump operator consists of the projector (1 — By)/2 
forming an interrogation part that checks whether the 
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FIG. 4: Sublattice schemes for the toric code. For the 
z = 2 partitioning scheme, the dominant contribution to the 
crosstalk occurs at a distance of v^a, while for the z — ^ 
scheme, the distance is increased to \/vfa. 

system is in a +1 eigenstate of By and a pump part ax\ 
which flips a single spin and thus turns any state with 
By = —1 into a state with By = +1. From the toric code 
Hamiltonian it is obvious that the jump operator chang- 
ing By = —1 to +1 will lower the energy of the system 
and thus constitute a cooling effect. Indeed, by perform- 
ing this cooling operation many times on all vertices and 
plaquettes, we can prepare the ground state of the toric 
code 

B. Full lattice model 

As shown above, it is possible to construct a quan- 
tum simulation of the toric code for both coherent and 
dissipative dynamics by iterating over all plaquettes and 
vertices in the full lattice system. For maximum per- 
formance, it is desirable to exploit the largest degree of 
parallelism, while keeping errors at a minimum. While it 
is possible to address many molecules in parallel at the 
same time, one has to be aware of a few limitations. First, 
one has to make sure that any parallel implementation 
of plaquette and vertex operators does not try to address 
the same molecule twice during a single operation. This 
is naturally enforced by operating plaquette and vertex 
operations sequentially and partitioning the system into 
two different sublattices = 2), see Fig. HI Second, the 
long-range tail of the dipole-dipole interaction will intro- 
duce errors due to crosstalk between molecules on dif- 
ferent plaquettes or vertices. This effect can be reduced 
by increasing the unit cell of the sublattice. Already by 
using a z = 5 sublattice, the crosstalk can be reduced to 
a percent level perturbation. 

C. Imperfections 

In any realistic experimental situation, there will be 
other sources of imperfections besides crosstalk between 
distant molecules due to the long range dipolar interac- 
tion, such as phase noise of the laser fields used for optical 
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FIG. 5: Dependence of the residual temperature T on the 
two qubit gate error probabihty e. Numerical simulation for 
a system of 16 ensemble molecules averaged over 1000 real- 
izations are shown as crosses, while the solid line is derived 
within linear response theory. The inset shows the relaxation 
dynamics for e — 0.01. 



addressing of individual molecules. Here, the dominant 
source of error will be related to the two qubit controlled 
phase gates. Effectively, such errors can be described in 
terms of a random phase flip acting on one of the 
qubits with probability e. These gate errors will lead to 
additional noise terms in the quantum master equation 
[l3|, and for the toric code will result in a finite anyon 
density n, which can be cast as an effective temperature 
T, according to 



2^0 



ks log n 



(9) 



Applying linear response theory to the quantum mas- 
ter equation ([8]), we can obtain an asymptotic expression 
for the anyon density n = 4c£, with the numerical con- 
stant c depending on the details of the error model. In 
the case of phase errors, we find c = 7/2. As shown 
in Fig. O this behavior is reproduced by numerical sim- 
ulations in the limit of small gate errors. Remarkably, 
already for gate errors on the order of £ ~ 0.01, the 
stationary state reaches a low temperature regime char- 
acterized by kBT/Eo < 1, where topological order can 
be detected within finite-size systems |40|. 



IV. EXPERIMENTAL REQUIREMENTS 

To be specific, we focus on an implementation using 
NaK molecules [4l|, which combine chemical stability 
(42I with a relatively large electric dipole moment of 
d = 2.7 D. Single qubit operations can be carried out 
with a speed of ~ 20 kHz and hence will be much 
faster than the two-qubit conditional phase gates. For 
the latter, we obtain a gate speed of ~ 2 kHz for 
an a = 532 nm optical lattice, according to Eq. (j3j). 
The overall timescale £^0 for the simulation of the toric 
code will be lower due to the subsequent implementa- 
tion of plaquette and vertex operators and the need for 
partitioning the system into several sublattices. Here, 
for the aforementioned z = 5 partition scheme, we ob- 
tain Eq = h X 40 Hz, which is compatible with typical 
timescales in experiments with ultracold polar molecules. 
For comparison, although the Rydberg quantum simula- 
tor can reach energy scales up to three orders of magni- 
tudes larger [l3|, the polar molecule quantum simulator 
can achieve similar results for coherence times on the or- 
der of tc ~ 100 ms, corresponding to a two qubit gate 
error of e ^ t^/tc = 0.005. 



V. CONCLUSION AND OUTLOOK 

In summary, we have demonstrated the feasibility of a 
quantum simulator for many-body spin interactions us- 
ing ultracold polar molecules trapped in optical lattices. 
The proposed architecture includes both coherent and 
dissipative dynamics and is quite robust against experi- 
mental imperfections, as we have exemplified for the toric 
code Hamiltonian. At the same time, the setup naturally 
allows for extensions to even larger classes of strongly cor- 
related spin models, paving the way towards the realiza- 
tion of a universal quantum simulator based on ultracold 
polar molecules. 
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